Beyond LISA: Exploring Future Gravitational Wave Missions 
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The Advanced Laser Interferometer Antenna (ALIA) and the Big Bang Observer (BBO) have been 
proposed as follow on missions to the Laser Interferometer Space Antenna (LISA). Here we study the 
capabilities of these observatories, and how they relate to the science goals of the missions. We find 
that the Advanced Laser Interferometer Antenna in Stereo (ALIAS), our proposed extension to the 
ALIA mission, will go considerably further toward meeting ALIA's main scientific goal of studying 
intermediate mass black holes. We also compare the capabilities of LISA to a related extension of 
the LISA mission, the Laser Interferometer Space Antenna in Stereo (LISAS). Additionally, we find 
that the initial deployment phase of the BBO would be sufficient to address the BBO's key scientific 
goal of detecting the Gravitational Wave Background, while still providing detailed information 
about foreground sources. 



I. INTRODUCTION 

The launch of the Laser Interferometer Space Antenna 
(LISA)[l| in the next decade will usher in a new era for 
gravitational wave detection. LISA will not only detect 
the presence of gravitational waves, but also provide de- 
tailed information about many thousands of gravitational 
wave sources, vastly expanding the burgeoning field of 
gravitational wave astronomy. 

The primary sources for LISA are expected to be com- 
pact galactic binaries 0, supermassive black hole bina- 
ries, and extreme mass ratio inspirals. The gravita- 
tional wave signals detected by LISA encode informa- 
tion about the parameters of the binary systems, such 
as sky locations, masses, and distances, which can be 
extracted by matched filtering. Studies into the accu- 
racy with which the source parameters can be recov- 
ered from the LISA data have been performed for com- 
pact galactic binaries [3, 0, IE 0- supermassive black 
hole binaries^ 0, SIS 13) an d extreme mass ratio 
inspirals[Ioj. 

It is hoped that LISA will be the first of several efforts 
to explore the low frequency portion of the gravitational 
wave spectrum accessible to space borne interferome- 
ters. One possible follow on mission is tentatively named 
the Advanced Laser Interferometer Antenna (ALIA) [ijj , 
featuring a spacecraft configuration similar to that of 
LISA, with smaller arm lengths and lower noise levels. 
The primary mission of ALIA will be detecting inter- 
mediate mass black holes (IMBHs), with masses in the 
range 50 — 50, 000Af Q . Information on populations, loca- 
tions, and event rates could greatly enhance theories of 
black hole formation and evolution. Anoth er p roposed 
follow on mission is the Big Bang Observer |12j (BBO). 
The BBO will be an extremely sensitive antenna that is 
designed to detect the Gravitational Wave Background 
(GWB) left by the Big Bang. According to the standard 
cosmological picture, the GWB is a relic of the early in- 
flationary period of the Universe. Just as the COsmic 
Background Explorer (COBE) and the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) missions provided 
information about the Universe around the time of last 



FIG. 1: The proposed orbital configuration of the Big Bang 
Observer. 



scattering, the BBO should be able to provide informa- 
tion about the earliest moments in the history of the Uni- 
verse. The current BBO proposal calls for four LISA-like 
spacecraft constellations, with the orbital configuration 
shown in Figure ^ Two of the constellations will be 
centered on a 20° Earth-trailing orbit, rotated 60° with 
respect to each other in the plane of the constellations. 
These constellations will be referred to as the star con- 
stellations, as the legs of the constellation sketch out a 
six pointed star. The remaining two constellations are to 
be placed in an Earth-like orbit 120° ahead and behind 
the star constellations. These two constellations will be 
referred to as the outrigger constellations. The purpose 
of the outrigger constellations is to provide greater angu- 
lar resolution for foreground sources (see Section 11111 for 
details). The star constellations provide maximum cross- 
correlation of gravitational wave signals between the two 
constellations, with minimal correlated noise [13]. The 
noise in each detector is expected to be independent, so 
over time the overlap of the noise between the two con- 
stellations will tend to average to zero, while the overlap 
of the signal will grow. The plan is to deploy the BBO in 
stages, starting with the star constellation, then adding 
the outrigger constellations at a later date. 

Figure |21 shows the detector sensitivities for LISA, 
ALIA, and BBO, a level for the expected extra-galactic 
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significantly improve the angular resolution for transient 
sources. This is followed by our study of the ALIA mis- 
sion in Section IIV Al and a study of our dual-detector 
variant of the ALIA mission in Section liVBI Results for 
unequal mass binaries for ALIA and ALIAS are in Sec- 
tion ^^21 The LISA and LISAS missions are compared 
in Section IIV Dl The BBO mission is studied in Sec- 



tion |IVEJ and results for a down-scoped version of the 
BBO mission are reported in Section TlV Fl Concluding 
remarks are made in Section IVI 



II. REVIEW OF PARAMETER ESTIMATION 



logl 0(f) (Hz) 

FIG. 2: Sensitivity curves for LISA, ALIA, and BBO with 
optimally filtered amplitude plots for equal mass binaries at 
z = 1 in their last year before coalescence. The masses shown 
are in units of solar masses. The squares denote the frequency 
one week from coalescence. 



background confusion noise 14], and optimally filtered 
amplitude plots for equal mass binaries in their last year 
before coalescence at redshift z = 1. The squares shown 
on the amplitude plots denote the frequency of the signal 
one week before coalescence. The signal from coalescing 
equal mass binaries at z = 1 with masses above 10 
will be detectable by both ALIA and BBO. The BBO 
will also be able to detect lower mass systems such as 
neutron star (NS) binaries. However, for coalescing bi- 
naries, a large portion of the signal strength is due to the 
rapid inspiral in the last week, and for 10M© binaries at 
z = 1 or greater, much of this power is deposited after 
the signal has crossed the ALIA sensitivity curve. 

Here we begin by analyzing the parameter resolution 
of ALIA for IMBHs. These studies indicate that by in- 
cluding a second identical constellation of spacecraft in a 
20° Earth- leading orbit, providing a constellation sepa- 
ration of 40°, the parameter resolution will be greatly 
improved. We will call this two constellation config- 
uration the Advanced Laser Interferometer Antenna in 
Stereo (ALIAS). We will also address the natural ques- 
tion, given the increases in precision that ALIAS pro- 
vides relative to ALIA, what would be gained by adding 
a second constellation to the LISA mission - the Laser 
Interferometer Space Antenna in Stereo (LISAS). Lastly, 
we analyze the parameter resolution of the BBO for NSs, 
10M Q , and 100M o black holes (BHs), and for the initial 
deployment phase of the BBO, which we call BBO Star, 
containing the two star constellations of the BBO con- 
figuration, but not the outrigger constellations. We will 
show that parameter resolution of the BBO and BBO 
Star should be sufficient to detect foreground binary sys- 
tems, out to and beyond a redshift of z = 3. 

This paper is organized as follows. We begin with a re- 
view of parameter estimation in Section^ In Section lTTTl 
we describe how a pair of well separated detectors can 



In our analysis, the orbits of the binaries are treated as 
quasi-circular and spin effects are neglected. This leaves 
nine parameters that will describe the binary systems: 
sky location (#, 4>); inclination and polarization angles 
(t,-0); reduced and chirp masses (/z, A4); time to coa- 
lescence (t c ), luminosity distance (-Dl), and the initial 
orbital phase (70). The signals are modeled using a trun- 
cated second-Post Newtonian (2PN) approximation fl5j 
whereby the amplitude is kept to Newtonian order while 
the phase is kept to second order. In other words, we only 
include the dominant second harmonic of the orbital fre- 
quency. 

The response of a space-borne instrument to a gravita- 
tional wave source is encoded in the Michclson-like time- 
delay interferometry (TDI) variables 01 > Xi(t). Here 
the subscript i denotes the vertex at which the signal is 
read out. In the equal-armlength limit, the TDI signal, 
Xi{t), can be formed from a time-delayed combination of 
Michelson signals, Mi(t), by 



X i (t)=M i (t)-M i (t-2L). 



(1) 



This differencing cancels the laser phase noise, while pre- 
serving the gravitational wave signal. Rather than work 
with the correlated Xi(t) variables directly, we use the 
orthogonal signal combinations |16| 



where 



A(t) 
E(t) 

Su(t) 



-(S n (t)-S n (t-L)), 



-(Sz(*)-5j(t-L)), 



{2M 1 (t) - M 2 (t) - M 3 {t)) , 



(2) 



-L(M 2 (t)-M 3 (t)) 



(3) 



The A, E combinations cancel the laser phase noise that 
would otherwise dominate the Michelson signals. At high 
frequencies, / > /* = c/(27ri), where L is the length 
of the detector arms, a third independent data channel 
T(t) becomes available. For frequencies below the T- 
channel is insensitive to gravitational waves and can be 
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FIG. 3: Sensitivity curves for the A, E, and T channels for 
ALIA with optimally filtered amplitude plots for equal mass 
binaries at z = 1 in their last year before coalescence. The 
masses shown are in units of solar masses. The squares denote 
the frequency one week from coalescence. 



TABLE I: SNRs for the T channel and the combined SNRs 
for the A and E channels as seen by ALIA for the equal mass 
binaries located at z = 1. 





51 Week SNR 
A+E T 


Last Week SNR 
A+E T 


1O 2 M 
1O 3 M 
1O 4 M 


44.6 0.599 
122 1.23 
131 3.52 


16.0 2.29 
302 16.1 
2145 34.8 



used as an instrument noise monitor. To simplify our 
analysis we did not use the T-channel. Including this 
channel would not have had much effect on the trend of 
our results. As an example, looking at Figure [3] one can 
see that for ALIA inclusion of the T-channel has little 
to no effect for binaries < 10 2 Af©. Higher mass bina- 
ries would have marginally improved signal strength and 
correspondingly improved resolution of the parameters. 
Table J] contains signal-to-noise ratio (SNR) information 
for the combined A and E channels, as well as the T 
channel for equal mass binaries as seen by ALIA in the 
year before coalescence. The SNRs in Table [TJ are for the 
final year before coalescence, separated into the contri- 
bution from the first 51 weeks and the contribution from 
the final week. 

In the low frequency limit, the variables Si and Su 
reduce to the independent data channels defined b y C ut- 
ler Q . We use the rigid adiabatic approximation |l9j ] to 
model the gravitational wave content of the data chan- 
nels. The instrument noise in the A and E channels was 
modeled by E3 



sin 2 (//2/»)[(2 + cos(//7,))S pos 

5accl 



S n (f) 

+2(3 + 2cos(///.) + cos(2///.)) 



(2tt/) 4 



/(2L) 2 



where S pos and S'acci are the one-way position and accel- 
eration noise contributions. Table ITU lists the instrument 
parameters used in our study of the LISA, ALIA and 
BBO missions. 



TABLE II: Instrument Parameters 



Parameter 


LISA 


ALIA 


BBO 


S pos (n^Hz- 1 ) 


4 x 10~ 22 


1 x 10~ 26 


2 x 10~ 34 


Saccoi (m 3 s- 4 Hz- 1 ) 


9 x 10" 30 


9 x 10" 32 


9 x 10~ 34 


L(m) 


5 x 10 9 


5 x 10 s 


5 x 10 7 



The A and E signals are functions of the nine param- 
eters A — > (6, <j>, t c ,DL,t, ^,A4,/x,7i) that describe the 
source. The source parameters can be recovered from 
the data channels using a variety of data analysis meth- 
ods. We will consider maximum likelihood estimation, 
using the log likelihood function 

log£(X') = (s\h(X'))-~(h(X')\h(X')), (4) 

where s denotes the data and h(X') denotes the search 
template. We use the standard noise-weighted inner 
product (a | b) summed over the independent data chan- 
nels: 



(a\b) = 2j2 



A,E 



a*(f)Hf)+a(f)b*(f) 

Sn(f) 



df. (5) 



For our work here the limits of the integration are set by 
the initial and cut-off frequencies of the binary systems. 
To first order these are given by: 



initial 



3.1 



/ 10 3 M Q \5/8 

\M(l + z)) 



ut-off 



4.4 



1(1 + z) 
10 3 M, 







mHz, 



Hz. 



(6) 



(7) 



^M tota i{l + z). 

where M total = .M 5 / 2 //// 2 is the total binary mass. Set- 
ting A' equal to the true source parameters yields the 
SNR 



SNR 2 = 2(log£(A)> = (h(X)\h(X)). 



(8) 



Here the angle brackets () denote an expectation value. 
The maximum likelihood estimator Aml is defined by: 



<91og£(A M L) 
dX i 



0. 



(9) 



The Fisher Information Matrix (FIM), T, is given as the 
negative of the expectation value of the Hessian evaluated 
at maximum likelihood: 



G> 2 log£(AM L ; 
dX l dXi 



(10) 



4 



where h t i = For large SNR, the parameter estima- 
tion uncertainties, AA Z , will have the Gaussian probabil- 
ity distribution 



detr 



1, 



p(AX) = J cxp J--r« AA*AA' ) , (LI ) 



and variance-covariance matrix 



(AA J AA J '> = C ij = (r- 1 ) 



(12) 



The uncertainties in each of the parameters are given 
by A A 1 = (C™) 1 / 2 (no summation). One can see from 
equations (J5J and (|12fl that the parameter uncertainties 
scale inversely with the SNR. 



III. MULTIPLE CONSTELLATION DETECTION 

Detectors such as LISA and ALIA determine the posi- 
tions of gravitational wave sources through both ampli- 
tude and frequency modulation. The angular resolution 
improves with time due to the accumulation of SNR and 
the synthesis of a long baseline as the detectors move in 
their orbit. In contrast, a detector array like the BBO has 
widely separated elements, and thus has a built-in base- 
line. Adding a second widely separated constellation to 
LISA or ALIA would increase the SNR by a factor of 
y/2, but the main gain in angular resolution for tran- 
sient sources would be due to the built in baseline. The 
advantage of a multi-element array can be understood 
from the following toy model. Suppose that we have a 
gravitational wave of known amplitude and frequency. 
Neglecting the effects of amplitude modulation we have 

h = Acos(2nf(t + Rsm(6) cos(2?r f m t + k- <f>))) (13) 

Here / is the source frequency, R is the distance from 
solar barycenter to the guiding center of the constella- 
tion, f m — 1/year is the modulation frequency, k is the 
azimuthal location (along Earth's orbit) of the guiding 
center, and 6 and (ft give the source location on the sky. 

With this two parameter (9 and <fi) signal we can an- 
alytically derive the uncertainty in the solid angle for a 
source observed for a time, T obs , by a single constellation 

25„(/) 



(AnfR) 2 sm(2e)T obs 
1 

x / 

v /(l-sinc 2 (2 7 r/ m T ofcs )) 



(14) 



For small observation times Af2 s i ng i e scales as T~ b 2 , while 
for large observation times it scales as T7* . 

Turning to the dual detector case, we simply add to- 
gether the FIMs for each individual detector. For two 
constellations separated by an angle Ak, this yields a 



solid angle uncertainty of 

Sn(f) 



Ail 



dual 



(AnfR) 2 sm(20)T obs 
1 



.(15) 



1 - sinc 2 (27r/ m T obs ) cos 2 (Ak)) 



For small observation times and non-zero An, Afiduai 
scales as T~j£. In other words, the built in baseline leads 
to a much improved angular resolution for short observa- 
tion times. This is very important for coalescing binaries 
as most of the SNR accumulates in the final days or weeks 
prior to merger. Note that if the two constellations are 
co-located, Ak = 0, the uncertainty in the solid angle is 
reduced by a factor of two relative to the single detector 
case by virtue of the increased SNR. Also, note that this 
toy model only includes Doppler modulation, thus the 
symmetry between An = and An = n. Including the 
amplitude modulation breaks this symmetry. 



IV. RESULTS 

The data shown here is for sources at a redshift of z = 1 
for ALIA, ALIAS, LISA, and LISAS and redshift z = 3 
for BBO, and BBO Star. These correspond to luminosity 
distances of 6.63 Gpc and 25.8 Gpc, respectively, using 
the best fit WMAP cosmology [2ij. Each binary system is 
observed for the last year before coalescence. Each data 
point is distilled from 10 5 random samples of 6, (ft, i, ip, 
and 7 . For t c , Dl, M., and \i we use logarithmic deriva- 
tives so that the uncertainties listed have been scaled by 
the value of the parameters. The uncertainty in sky lo- 
cation is simply the root of the solid angle uncertainty 
(V Ail). The remaining angular parameters, i, tft, and 
7o, have not been scaled. 

The angular variables were chosen by using a Monte 
Carlo method. The values for cos(0) and cos(t) were cho- 
sen from a random draw on [—1, 1]. Values for (ft, ip, and 
7o were each chosen from random draws on [0, 2tt\. The 
parameters A4 and Dl were set for each Monte Carlo 
run (though they changed between runs). Time to coa- 
lescence, t c , was set to 1 year plus a small offset so that 
during the year of observation the binary did not reach 
a relativistic regime that would not be properly modeled 
by the 2PN approximations used. 

The mean SNRs quoted in this paper are calculated 
by taking the square root of the average of the squares of 
the individual SNRs. For our analysis, positive detection 
will be restricted to SNRs above 5. 



A. Results for ALIA 

Table IIIII summarizes the medians and means of pa- 
rameter uncertainties for detections by ALIA of equal 
mass binaries with masses of 10 2 M P) , 10 3 M(7,, 10 4 M, 
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FIG. 4: Histograms of the parameter uncertainties and SNR, 



for equal mass binaries of 10 M® 
ALIA. 



at z = 1, as detected by 



and 10 5 Mq. The SNRs for this range of masses shows 
that ALIA should get positive detection of 99 + % of co- 
alescing IMBHs located at z = 1 or closer, which would 
provide good information on the coalescence rates. The 
great precision in the measurement of M. and /x will 
provide a clear picture of the constituent masses of the 
binary systems. Furthermore, the sub-degree precision 
in the sky location, combined with luminosity distances 
known to a few percent, will facilitate the construction 
of a three dimensional distribution of IMBHs with which 
to test theoretical predictions (see Ref. [2lJ for a related 
discussion concering LISA and supermassive black holes). 

The uncertainties in t c a month from coalescence (using 
11 months of data) will be on the order of a few minutes. 
This will provide warning time for ground-based gravi- 
tational wave detectors, as well as other systems (tele- 
scopes, neutrino detectors, etc.) to gather as much and 
as varied data as possible about the coalescence. 

The shapes of the histograms shown in Figure^ which 
are from the 10 3 M Q data, are representative of the his- 
tograms for each of the detectors covered in this work. 
Note that for l, ip, and 7 , the tails of the histograms 
run far beyond the range of the plots shown, raising the 
values of the means considerably above their respective 
median values. For example, while the median value of 
7o for equal mass binaries with masses 10 2 M Q is 12.2° 
its mean value is 2185°, which is well beyond the [0, 27r] 
range of 7 Q . Uncertainty ranges that are larger than the 
possible range of the parameter tell us that the param- 
eter is indeterminate. When uncertainty ranges exceed 
the possible range of a parameter, one may drop that 
parameter from the FIM analysis. For our analysis, we 
did not discard any parameters. As an example of how 
this affects the analysis, consider the 1O 2 M study where 
the mean value of the uncertainties that lie in the [0, 27r] 
range for -f a would be 36.3°, while its median would be 
9.79°. However, for 16.9% of the binaries, the 7„ param- 
eter would be indeterminate. 



Figure [5] shows the SNR histograms for 10M© equal 
mass binaries for ALIA (and ALIAS). As was expected 
from Figure [3 the SNR values for ALIA are low, with 
nearly 60% too low for a positive detection. Of note, 
though, is that the SNR scales roughly as the luminosity 
distance. This suggests that ALIA should be able to 
detect nearly all IOMq binaries with luminosity distances 
less than 2 Gpc. For binaries with masses beyond 15Af Q 
the low end of the range of SNRs is above 5, meaning that 
equal mass binaries with masses in the range of IMBHs 
should be detectable out to z = 1. 



B. Results for ALIAS 

The main purpose of ALIA is gathering information 
about IMBHs; our data shows it is capable of doing this 
with some success. A more accurate IMBH census could 
be derived from a dual constellation version of ALIA that 
we call ALIA in Stereo or ALIAS. Each component of the 
ALIAS constellation would be offset from the Earth by 
20°, one in an Earth-trailing and the other in an Earth- 
leading orbit, giving a 40° separation in order to provide 
increased parameter resolution in the IMBH range. 

Table llVl summarizes the medians and means of the pa- 
rameter resolutions that could be achieved by the ALIAS 
mission. Results are given for equal mass binaries with 
masses of 1O 2 M , 10 3 M Q , 10 4 M Q , and 10 5 M Q . As ex- 
pected, the mean SNR increases by \[2 relative to ALIA. 
The improvements in parameter resolution are, however, 
considerably larger. At the upper end of the IMBH mass 
range the angular resolution improves by a factor of ~ 90 
and the luminosity distance resolution improves by a fac- 
tor of - 23. 

For masses below ~ 10 2 M Q , ALIAS provides roughly 
a factor of \[2 increase in the median values of each of the 
parameter uncertainties. This can be seen in the 10 2 Af Q 
data in Tables IIIII and IIVI as well as in Figure [fj] and 
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FIG. 5: SNR histogram for ALIA and ALIAS from 10M© 
equal mass binaries at z — 1. 
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TABLE III: SNR and parameter uncertainties for ALIA with sources at z = 1. 





1O 2 M 


1O 3 M 


1O 4 M 


10 s M e 




median 


mean 


median 


mean 


median 


mean 


median 


mean 


SNR 


71.7 


76.2 


466 


567 


2937 


3641 


6481 


7970 


Sky Loc. 


0.169° 


0.201° 


0.234° 


0.264° 


0.382° 


0.417° 


0.518° 


0.546° 


ln(tc) x 10" 8 


1.38 


1.58 


2.74 


3.11 


4.41 


5.16 


5.38 


6.47 


ln(M) x 10 s 


9.56 


11.0 


12.1 


13.9 


31.3 


35.7 


136 


152 


ln(^) x 10" 5 


46.8 


53.6 


5.10 


5.86 


1.76 


2.03 


2.66 


3.03 


ln(D L ) x 10 -2 


6.25 


32.7 


1.34 


6.44 


1.08 


3.50 


1.15 


3.08 


i 


3.99° 


125° 


0.762° 


22.4° 


0.513° 


8.00° 


0.515° 


5.66° 


i> 


5.44° 


1092° 


1.08° 


176° 


0.840° 


42.3° 


0.931° 


24.7° 


Jo 


12.2° 


2185° 


4.15° 


352° 


2.79° 


85.4° 


2.51° 


50.0° 



Figure which plot the uncertainty in the sky location 
and luminosity distance, respectively, against the chirp 
mass of the binaries. The trend shown in these figures 
holds true for D^, i, and ip in that mass range, while 
the increases in t c , /z, j Q , and M. are slightly larger, 
~ 1.55 better than the values for ALIA. However, the 
median SNRs are below 5 for equal mass binaries with 
chirp masses below ~ 8M© for ALIAS (below ~ IOMq 
for ALIA). 

Figure [S] also shows that the benefits of a dual con- 
stellation becomes even more significant when the final 
chirp of the binary occurs above the detector noise (see 
Figure |2| . The main improvement is in the angular res- 
olution, but the decreased covariances between the sky 
location and other parameters lead to improved measure- 
ments of t c , Dl, t, and i[> at the upper end of the IMBH 
mass range. The resolution of these parameters improves 
(relative to ALIA) by a factors of 33, 23, 18, and 21, re- 
spectively. The parameters M., /i, and 7 Q only see a 
factor of ~ 2 improvement in resolution. 

Comparing Figure with Figure H3 one can see that 
the increased resolution in the luminosity distance corre- 
sponds with the increased precision in sky location due 
to their large covariance. A similar effect was seen in 
the work of Hughes and Holz [z| , where the addition of 
electromagnetic information to fix the sky location of a 
BH merger resulted in a marked increase in resolution of 
the luminosity distance. In our case, the degeneracy is 
broken by the improved angular resolution afforded by 
the baseline between the ALIAS detectors. These im- 
provements in sky location and luminosity distance res- 
olution will provide an even more detailed three dimen- 
sional distribution of the IMBHs than that provided by 
ALIA. Increasing the 3-D resolution of the distribution of 
IMBHs by roughly a factor of 6 on the low end of the 
mass range (50M Q ) up to ~ 200,000 on the high end of 
the mass range (50,000M Q ), would provide much more 
detailed information compared to ALIA, with which to 
test theoretical predictions about IMBHs. 

Figure shows the SNR histogram for 10Af© equal 
mass binaries for ALIA and ALIAS. While more than 
two-thirds of the sources at z = 1 would be positively 
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FIG. 6: Uncertainty in the sky location for both ALIA and 
ALIAS plotted against binary chirp mass. Binary system con- 
stituents in this plot have equal masses and are located at 

2 = 1. 
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FIG. 7: Uncertainty in the luminosity distance of both ALIA 
and ALIAS plotted against binary chirp mass. Binary system 
constituents in this plot have equal masses and are located at 

2 = 1. 
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TABLE IV: SNR and parameter uncertainties for ALIAS with sources at at z = 1. 
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1.34° 
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detected, nearly one-third would be missed. However, 
these results suggest that ALIAS will be able to detect 
nearly all of the IOMq binaries with luminosity distances 
less than 3 Gpc. This is roughly in keeping with the 
general trend of ALIAS'S capabilities compared to ALIA. 

From Figure [21 we saw that for masses below IOMq, 
the source's signal lies below the sensitivity curve for the 
week prior to coalescence. For larger masses, more of 
the binary's final chirp will be detectable. This is why 
the angular resolution of ALIAS becomes so much better 
than ALIA's at the upper end of the IMBH mass range. 
For masses above 10 5 M Q the final chirp occurs before the 
"sweet spot" of the sensitivity curve, which diminishes 
ALIAS'S advantage over ALIA. 



C. Results for ALIA & ALIAS for unequal mass 
binaries 

Equal mass binaries are the easiest to study, but they 
can give an overly optimistic picture of the instrument 
capabilities as they yield the smallest parameter uncer- 
tainties. To study this bias and provide a more realistic 
picture of the capabilities of ALIA and ALIAS, we now 
consider the case where the masses of each component in 
the binary are drawn from logarithmic distributions in 
the range 1 — 10 8 M Q . The results are still presented as 
a function of chirp mass, but the mass ratios reflect the 
underlying mass distribution. 

Figure 0D shows the plot of sky location uncertainty 
against chirp mass. The uncertainties in sky location 
for masses above 10 2 M Q are larger than those shown in 
Figure EJ The increase in sky angle uncertainty relative 
to the equal mass case was a factor of ~ 2 for ALIA. 
Similarly, other parameters show a factor of ~ 3 increase 
in their uncertainties relative to the equal mass study. 

However, ALIAS still maintains an advantage in locat- 
ing binaries. In fact, for chirp masses below ~ 500M Q the 
increase in precision for ALIAS over ALIA in locating un- 
equal binaries is nearly the same as it was for equal mass 
binaries, a factor of roughly ^/2. At higher chirp masses 
we see less of a difference in angular resolution than in 
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FIG. 8: Uncertainty in the sky location for both ALIA and 
ALIAS plotted against binary chirp mass. Binary system con- 
stituents in this plot have unequal masses and are located at 

2 = 1. 



the equal mass study (factors of ~ 9 compared to ~ 90). 
This trend holds for t c , Dl, i, and ip, which show max- 
imum increases in precision ~ 4 for ALIAS over ALIA, 
while M } n, and 7 D show maximum increase in precision 
~ 2 for ALIAS over ALIA (as they did for equal mass 
binaries). Thus the increased resolution in the 3D distri- 
bution of IMBHs for ALIAS over ALIA for unequal mass 
binaries ranges from a factor of ~ 6 up to ~ 300. 



D. Results for LISA and LISAS 

Figure [5] compares the sky location uncertainties for 
the LISA mission to the LISAS mission. As was seen 
with ALIA and ALIAS, the addition of a second con- 
stellation to LISA provides a marked increase in param- 
eter resolution. The increases in parameter resolution 
for LISAS over LISA are similar to those found between 
ALIAS and ALIA. The angular resolution showed a max- 
imum improvement of ~ 25 just above 10 5 Mq. Also, a 
lesser benefit than that shown in Figure |H1 occurs in the 
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FIG. 9: Uncertainty in the sky location for both LISA and 
LISAS plotted against binary chirp mass. Binary system con- 
stituents in this plot have equal masses and are located at 
z = l. 
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FIG. 10: Uncertainty in the luminosity distance and the SNR 
for both LISA and LISAS plotted against binary chirp mass. 
Binary system constituents in this plot have equal masses and 
are located at z — 1. 



t c , Dl, i, and ip parameters (with factors of maximum 
increase ~ 12 just above 1O 5 M ), while M., fi, and j 
show only a modest maximum improvement by a factor 
of ~ 2. Our results are in the same range as those found 
by Seto |2^| for sources detected by LISA that have un- 
dergone strong gravitational lensing. In that case the 
extended baseline was provided by the time delay in the 
arrival of the signals, which effectively turned a single 
LISA detector into a LISAS system. Seto also gave brief 
consideration to the performance of a dual LISA mission, 
and found the ratio of angular resolution between LISAS 
and LISA to be larger than that seen in our simulations. 
The discrepancy can be traced to Seto including a large 
galactic confusion background which limits the time that 
the sources are in-band, thus magnifying the advantage 
of the dual configuration. 

The shift in the placement of the minimum uncertain- 



ties, compared to ALIA and ALIAS, can be understood 
from Figure [3 The signal from a 10 3 M Q equal mass bi- 
nary is passing below the LISA sensitivity curve in the 
last week before coalescence, and before the final chirp of 
the binary. This is evident in Figure^] where the angular 
resolution of the single and dual configurations is similar 
for masses below 10 3 M©. The sweet-spot of the LISA 
sensitivity curve is at ~ 5 mHz, which is where 10 Mq 
binaries experience their final chirp. Above 10 5 M Q , less 
of the final chirp is occurs near the sweet spot, and the 
difference in the angular resolution diminishes. Figure lTUI 
shows how this shift corresponds precisely to the shift 
in the maximum SNR of the LISA and LISAS missions. 
The plot shows how the increase in the mean SNR (V2) is 
uniform across the range of chirp masses. Also, it shows 
again how the resolution in Dl improves with the addi- 
tion of more information in the form of increased angular 
resolution. 



E. Results for BBO 

Table summarizes the medians and means of the 
parameter uncertainties for detections of equal mass bi- 
naries with masses of 1.4M , IOMq, and 1O 2 M for the 
BBO. Histograms for this data have the same general 
shapes as those shown in Figure0]and so are not shown. 

As can be seen from the data, the BBO is an extremely 
sensitive detector. The SNRs show that if the BBO meets 
design specifications there will be positive detection of co- 
alescing stellar mass binaries out to (indeed well beyond) 
a redshift of z = 3. Similarly, the BBO will provide a 
precise picture of coalescence rates, including how these 
rates relate to luminosity distance. With a combination 
of sky location and luminosity distance the BBO should 
be able to pick out the host galaxy of the majority of 
coalescing binaries. Data taken by the BBO up to a 
month before coalescence (using 11 months of data) will 
determine t c to within seconds, again providing ample 
warning time for other detectors to gather data on the 
coalescence. 



F. Results for BBO Star 

As can be seen in Figure |2 the final chirp of solar 
mass or NS coalescing binaries occurs below the noise 
level of the BBO. Thus the outrigger constellations do 
not provide an extended baseline for the last few days of 
NS chirp. Positive detection of NS binaries at z — 3 can 
be accomplished with the initial deployment of the two 
constellations that make up the star, while still providing 
the cross-correlating needed to detect the GWB. 

Table IVII summarizes the medians and means of the 
parameter uncertainties for the detection of equal mass 
binaries with masses of 1.4M , IOMq, and 1O 2 M for 
the BBO Star. Histograms for this data have the same 
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TABLE V: SNR and parameter uncertainties for BBO with sources at z = 3. 
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median 
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ln(t c ) x 10~ 12 
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200 
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5.70 
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28.9 
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57.3° 
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96.1° 


0.0806° 


15.2° 


7° 


7.71° 


996° 


3.93° 


194° 


2.82° 


32.3° 



general shapes as those shown in Figure 0] and so are not 
shown. 

Figure ITTI plots the uncertainty in sky location against 
chirp mass for the BBO and BBO Star. As can be seen, 
the effect of the outrigger constellations is significant, 
providing for increased parameter resolution for the full 
BBO throughout the range of chirp masses shown, with 
a maximum increase of ~ 3700 around 200-Mo. Similar 
to the ALIA/ ALIAS comparison this extra precision in- 
creases until all of the final chirp lies above the sensitivity 
curve (see Figure |5Jl, and begins to decrease as less and 
less of the final chirp occurs in the frequency range of the 
sweet-spot of the sensitivity curve (which for BBO is ~ 1 
mHz). The improvement in the mean SNR for BBO over 
BBO Star is the expected y2- Also, the improvement in 
the resolution of t c (~ 3550) is comparable to that seen 
in the angular resolution. More modest increases occur 
for Dl, i, and ip, which increase by factors of 52, 35, and 
47, respectively. Only slight (~ 2.5)increases are seen for 
M, (J,, and j . 

The SNRs from BBO Star are sufficient for positive 
detection of binaries with constituents less than one so- 
lar mass out to, and beyond, z — 3. While the full BBO 
offers considerable adavantages for doing precision gravi- 
tational wave astronomy, BBO Star could fulfill the main 
science objective of detecting the cosmic gravitational 
wave background while still providing useful information 
about binary populations. 

V. CONCLUSION/DISCUSSION 



return a far more accurate census of the IMBH popula- 
tion. In addition we have shown that a similar extension 
to LISA would greatly improve its ability to locate the 
host galaxies of coalescing binaries. On the other hand, 
if we are willing to give up some of the precision astron- 
omy offered by the full BBO, we could get by with just 




Chirp Mass (in M sun ) 



FIG. 11: Uncertainty in sky location versus chirp mass for the 
standard BBO design and the star constellations of the BBO. 
Binary system constituents in this plot have equal masses and 
are located at z = 3. 

the first phase of the BBO deployment. The BBO Star 
configuration could satisfy the primary goal of detecting 
the CGB, while still providing a detailed binary census. 



While our survey is by no means comprehensive, it has 
helped to map out the science that can be done with the 
ALIA and BBO missions. We have shown that ALIAS, a 
modest extension to the ALIA mission, would be able to 
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